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ANALYTICAL  DESIGN  FOR  INTERNAL  BURNING  STAR  GRAINS  OF  SOLID  ROCKETS 


Lii  Chang-tang 

With  regard  to  the  design  of  internal  burning  star  grains  of 
solid  rocket  engines,  they  have  been  frequently  designed  based  on 
pure  geometrical  studies  and  their  corresponding  trial  methods.  The 
formula  and  computation  curves  shown  in  [1]  are  representative  which 
are  widely  used  in  engineering  design  and  scientific  studies  both  in 
our  country  and  abroad.  Actual  practice  showed  that  the  pure  geometri 
cal  formula  and  curves  given  in  [1]  are  not  only  too  complex,  but  also 
incomplete.  The  corresponding  trial  method  not  only  involves  a  huge 
computational  load,  but  also  is  rather  blind.  It  is  very  difficult 
to  ensure  various  technical  objectives.  Especially,  it  is  extremely 
difficult  to  reach  the  optimum  design.  In  this  paper,  we  will  attempt 
to  use  analytical  design  to  replace  trial  design.  Different  from  the 
traditional  pure  geometrical  studies,  we  combined  the  various  geometri 
cal  parameters  of  the  internal  star,  the  various  characteristic  para¬ 
meters  of  the  solved  propellent  and  the  technical  objectives  of  the 
engine  to  establish  the  grain  design  equation  series  according  to  the 
optimization  principle  of  grain  design.  It  was  matched  with  simple 

t 

computational  formulas  and  curves  to  quickly  solve  the  equation  series 
Compared  to  [1-3],  not  only  the  computational  load  is  less,  but  also 
the  required  technical  objectives  are  assured.  In  addition,  it 
approaches  the  optimum  design  in  improving  the  performance  of  the 
engine  so  that  significant  savings  in  materials  and  ease  of  computer 
calculation  are  realized. 

9  1 .  Geometrical  functional  relationship  of  the  Internal  star.  As 

shown  in  Figure  1,  based  on  the  Plobert  combustion  law,  it  is  not  dlff 
icult  to  obtain  the  functional  relations  between  the  geometrical 
parameters  of  internal  star  grains  presented  in  Figure  1.  To  simplify 
the  calculation*  only  functions  k  (Figure  2)  and  s^  (the  definitions 
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are  shown  in  Figure  1),  which  are  related 
to  the  star  corner  number  N  and  the  half 
angle  of  the  star  tip  0Q/2  are  given: 


"neN 

*•  —  2 A7  -  — 2-  +  2*/(l  —  e)  +  */ 
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Therefore,  the  first  stage  combustion 
Figure  1.  Diagram  of  area  Afel  and  the  initial  gas  passage  area 

inte?nainsta?f  ^  Apol  can  be  slmPlified  into  the  tw0  following 

equations: 

L,U  +  V.  +  2*0-/,)]  (3) 

—  tj  +  Pc  +  *fi  —  o  Wf  +  /{)  ^  ^ 

where  Lp  is  the  length  of  the  grain  and  other  symbols  are  as  presented 
in  Figure  1.  The  value  of  c  is  given  in  the  following  equation: 

*"(!"“  ^)»  +  *V  sin  (  ro>(  —  .V  fin1*  —  etc  —  ^  5  ) 
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Based  on  analysis,  we  know  that  the  angle  0/2  varies  in  a  mono- 
tonically  increasing  manner  from  the  star  side  disappearing  point 
y0Q/2  with  combustion  time  in  the  region  of  [0q/2,  0^2].  It  is 
called  the  dynamic  angular  variable.  0w/2  is  the  star  tip  half  angle 
when  combustion  stops.  0/2  has  the  following  relation  with  y  which 
varies  with  time: 


<j  +  /)// 


Using  this  dynamic  angular  variable  0/2,  the  function  of  the 
second  stage  combustion  area  A^H  and  residual  grain  area  A f  are  sim¬ 
plified  into  the  following  forma: 


izw+vr 
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(8) 
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where  the  residual  grain  area  star  tip  half  angle  0^2  when  combustion 
stops  is 
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Prom  mathematical  analysis,  we  know  tbat.A^^  has  a  minimum  Ahm^n 
at  5/2  and  its  dimensionless  form  is  (Figure  3  is  the  dimensionless 
functional  curve) 


2  N 
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Similarly  it  can  be  proven  that  Af  has  a  minimum  Afinin  at 
-i-.i,  i.e., 
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The  existence  of  *\min  and  Af.min  is  an  important  characteristic 
of  internal  star  grain  cylinders. 


The  effective  grain  area  kQQ  without  star  tip  angle,  i.e.,  f^O, 
is  very  important  to  the  establishment  and  solution  of  the  design 
equation.  It  is  given  by  the  following  equation: 

+  NPa$,&d*,$  —  —  N(ir  +  «i»~rco,'T') 
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Therefore,  when  there  are  star  tip  angles,  i.e.,  f^O,  the  effective 
grain  area  A@ol  is 

-J^r—A*-  «.H2*  -  OR  (13) 

The  minimum,  radius  RQ  of  the  internal  star  at  f^*0  provides  the 
constraint  relation  between  the  star  shaped  geometrical  parameters. 
Furthermore,  because  it  la  always  true  that  ALMJt'lrr  ,  R  is  a 


decreasing  function  of  the  angular  coefficient  e.  Let  RQ>0,  then 
the  upper  limit  of  e  which  is  en  upper  is: 

eN  upper  —  H  (aresin  t  •+  ( 
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As  long  as  e<eXI  .  .  then  R  >0.  Furthermore,  R  ,  >0  too  (when 

°  N  upper  o  ol 

fj/0).  . 

2.  The  establishment  of  the  grain  design  equation  series 

The  known  conditions  are: 

1)  the  technical  objectives  required  by  the  engine:  the  total 
impulse  value  Ir  or  in  the  required  range,  the  thrust  law  pmin”pmax> 
the  working  time  tmin-t  ,  the  working  temperature  TTnin"‘Trnay0c»  the 
external  diameter  of  the  engine  D  and  other  limitations. 

2)  the  characteristics  of  the  solid  propellent:  the  specific 
impulse  Isp,  the  combustion  speed  r  — the  temperature  coefficient 

the  flow  sensitivity  coefficient  limit  K^,  the  thrust  coefficient 
the  characteristic  speed  C  ,  the  critical  pressure  P^p,  and  the 
density  p. 

The  design  equation  series  obeys  the  grain  design  optimization 
principle:  to  ensure  the  realization  of  the  total  impulse  and  thrust 
requirements  as  well  as  the  working  time,  to  carry  out  stable  and 
normal  combustion,  to  have  a  large  effective  grain  coefficient  ne  and 
£  small  residual  grain  coefficient  nf.  Stress  is  concentrated  in  the 
allowable  range . 

In  order  to  establish  the  design  functions,  we  must  carry  out 
the  following  transformation: 

Transform  the  total  lmpuse  into  the  volume  of  the  affective 
grain  quantity 


184 
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The  thrusts  are  transformed  into  grain  combustion  areas 


c*p  -  cP'  - 
e*p  \ePm  /+: 


where  At  is  the  critical  area  of  the  nozzle,  a,  —  F.i./p„i.c,_rc, with >  /*«p. 

The  working  time  of  the  engine  is  transformed  into  the  combus¬ 
tion  arch  thickness  W  which  is  calculated  from  ir  —  rt  or»'  —  2V,/(.4kmim  + 
and  also  verified  experimentally. 

The  grain  cylinder  diameter  Dj.  is  D  minus  the  thicknesses  of  the 
shell,  gap  and  packaging. 

After  W  is  determined,  in  order  to  make  Af  approach  Af  min,  the 
star  angle  transition  arc  radius  f  is  calculated  according  to  the 
reference  value  f  and  the  adjustment  formula  ctf. 


U  —  0.25D,  —  w 

(15) 

„ _ w  -  fc) 

0.250,  -  (/  -  /,) 

(16) 

Adjust  f  to  make  |of|  small  to  the  extent  possible.  Simultaneous¬ 
ly,  under  conditions  allowable  by  grain  stress,  small  f  should  be 
abosen  in  order  to  obtain  better  design  performance. 

Therefore,  the  design  equation  to  ensure  the  engine  total  impulse 

requirement  is  _ _ _ 

v.  (17) 

The  design  function  t»  ensure  the  lower  limit  of  thrust  is 
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The  initial  pressure  peak  limit  equation  to  ensure  stable  combustion 


A+  +  0.5C2»  -  0/5 


(19) 


Eliminate  the  grain  cylinder  length  Lp  and  f1  from  the  above 
series  of  equations  and  then  substitute  equations  (1),  (2),  (8),  (12) 
into  them.  We  can  obtain  the  equation  related  to  the  angular  coeffi¬ 
cient  e  as: 
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(20) 


We  can  see  that  this  equation  simultaneously  includes  the  re¬ 
quired  technical  objectives,  the  propellent  characteristics  and  the 
internal  star  geometrical  parameters. 

3.  The  solution  of  the  design  equation  and  the  acquisition  of 
the  grain  geometrical  parameters.  Obviously,  it  is  difficult  to 
directly  solve  for  e  from  equation  (20).  For  this  purpose,  let  us 
change  e  into  two  functional  curves  Y1  and  Yg. 


where 


Y,  -  (i  DJ  -  A,  -  y,  -  L,u  -  “>/ 2A(2w  -  O  1 

A  *  Am  -  WL,) 


(21) 

(22) 


A  is  the  necessary  formula  to  Judge  the  solution.  Only  where 
A>0  there  is  a  solution. 

• 

For  each  positive  integer  N,  there  is  an  upper  limit  0N  Upper^2 
for  the  angle  0o/2  which  makes  A*0  and  f^-0.  0N  Upper/2  is  an  increas¬ 
ing  function  of  N;  therefore  0o/2  must  he  determined  in  the  range  of 
(0*  0jj  uppap/2)  in  combination  with,  the  initial  thrust  requirement. 

N  and  9-/2  are  the  variable  parameters  of  equation  (20).  In  order  to 
© 

ham  bet  tear  overall  performance  for  the  grains,.  K  and  0Q/2  should  he 
chosen  to  he  sm»  1  (e.g.  0o/2»lO° ) .  After  N  and  0Q/2  are 
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Figure  2.  The  k  function  curves  used  in  calculation 


Figure  3»  The  curvee  of  the  function  used  in 

calculation  mj'n 


selected,  use  (N^0Q/2)  to  check  the  k  value  In  Figure  2,  and  then 
calculate  the  functional  curves  Y-^  and  Y2  of  e  .  For  each  chosen  e , 
use  (N,e)  to  look  up  the  Sm^n/1  value  in  Figure  3.  Substitute  into 
the  following  equation  to  obtain  L  : 

r*' 


Consequently,  V  /L  can  be  obtained. 

e  p 

into  equation  (12)  to  get  A  in  order  to  obtain  A. 

eo 


(23) 

Substitute  e  and  N,  9Q/2  and  k 


equations  (2)  and  (8),  we  can  obtain  sQ  and  Ap. 


When  A>0 ,  from 
Furthermore ,  from 
corresponding  to  e . 


equation  (21)  we  can  get  the  functions  Y^  and  Y^ 

Similarly,  a  series  Y^  and  Y2  functions  of  e  are  calculated.  The 
abscissa  of  the  intersect  point  of  the  decreasing  function  Y^  and  in¬ 
creasing  function  Y2  of  e  gives  the  solution  es0iut-£0n  of  equation  (20 
Therefore,  the  angular  coefficient  of  the  grain  is  obtained.  The 
ordinate  is  the  ApQ1  of  the  grain.  From  (N,  eSO]_ution) »  we  can  obtain 
S^in/1  by  checking  Figure  3*  Similarly,  to  the  procedures  described 

\T  /  T  A  A  A 

f  * 

culated. 


before,  the  grain  length  Lp,  Ve/Lp,  Ago,  A,  Af,  sQ,  etc.,  can  be  cal- 


V  /L  is  the  value  of  Aa  , . 
e  p  eol 


Use  the  following  two  equations  to  calculate  fx  and  RQ1 


(24) 

(25) 


When  fi  =  0,  from  equation  (14)  we  know  that,  as  long  as  e<en  upp( 
then  it  is  necessary  to  have  Ro>0.  When  f^O,  as  long  as  the  equation 
has  a  solution,  Ro^>0  is  valid.  Furthermore,  the  smaller  0q/2  is, 
the  larger  Rq1  becomes. 


Discussion  on  the  solution:  When  varying  e  and  A<0  still  exists, 
i.e.,  the  equation  has  no  solution,  it  is  necessary  to  choose  a  larger 
N  and  a  smaller  0q/2.  It  is  also  possible  to  vary  Kkp  and  other  known 

parameters  in  the  allowable  range  (Note:  larger  N  means  larger  Rq1> 
larger  02/2  means  smaller  A^.  When  necessary,  we  must  adjust  N  and 

0q/2  simultaneously) , 


Up  to  this  point,  the  engire  geometrical  parameters  of  Internal 
star  grains  are  solved*  Finally,  let  us  calculate  the  variation  of 
the  combustion  surfaces 
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CO,  y0  /2]  is  the  first  stage  of  combustion  Inntjj 

which  is  the  linear  stage,  y — - z - / 

CO*— 3 
.2 

In  the  interval  tO,f,],  we  have 

A* I -M/|  +  2*0 '"/*)!*■*  (26) 

In  the  interval  [f^,  y0Q/2],  we  have 

A„  -  0.  +  ky)L,  (27) 


l yf/i»  w  1  is  the  secon<*  stage  of  combustion  which  is  the  non¬ 
linear  combustion  stage.  Using  the  dynamic  angular  variable,  0/2,  we 
can  quickly  obtain  the  variation  of  Afell  from  equation  (7). 


Incidentally,  P.  A.  Williams  et  al  stated  in  [2]  that  "in  the 
second  st-4ge  of  combust  ion. ..  the  grains  are  the  increasing  type,  i.e., 
increases  linearly  according  to  y=rt".  This  is  not  accurate. 

In  the  second  stage,  Abll  is  nonlinear;  furthermore,  when  0Q/2  <  5/2, 
it  decreases  before  it  begins  to  increase. 
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